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CN ■ Abstract 

The problem 2-LOCAL HAMILTONIAN has been shown to be complete for the quantum 
computational class QMA [TJ. In this paper we show that this important problem remains 
QMA-complete when the interactions of the 2-local Hamiltonian are between qubits on a two- 
dimensional (2-D) square lattice. Our results are partially derived with novel perturbation 
\Q • gadgets that employ mediator qubits which allow us to manipulate fc-local interactions. As a side 

^ \ result, we obtain that quantum adiabatic computation using 2-local interactions restricted to a 

2-D square lattice is equivalent to the circuit model of quantum computation. Our perturbation 
method also shows how any stabilizer space associated with a fc-local stabilizer (for constant fc) 
can be generated as an approximate ground-space of a 2-local Hamiltonian. 
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*0 ' 1 Introduction 

o . 

The novel possibilities that quantum mechanics brings to information processing have been the 
Q_i| subject of intense study in recent years. In particular, much interest has been devoted to under- 
standing the strengths and weaknesses of quantum computing as it pertains to important problems 
in computer science and physics. 

An important part of this research program consists of understanding which families of quantum 
systems are computationally complex. This complexity can manifest itself in two ways. On the one 
hand, a positive result shows that a given family of systems is "complicated enough" to efficiently 
^ . implement universal quantum computation. On the other hand, a negative result shows that 
certain questions about such systems are unlikely to be efficiently answerable. A proof of QMA- 
completeness offers compelling evidence of the negative kind while also locating the given problem 
in the complexity hierarchy, since QMA, -the class of decision problems that can be efficiently 
solved on a quantum computer with access to a quantum witness-, is analogous to the classical 
complexity classes NP and MA. More precisely, the class QMA is defined as 

Definition 1 (QMA) A promise problem L = L yes U L no C {0, 1}* is in QMA if there is an 
efficient (of poly(|a;|) size) uniform quantum circuit family {V x } xe {Q iy such that 

Vx e L yes , 3 IVfe) 6 H® pol ^\ Prob(F x .(|^)(Vx|) = 1) > 2/3, 

and 

VxeL no , V|0 eW 0poly(N) , Prob(^(|0(£|) = 1) < 1/3. 



*IBM Watson Research Center, Yorktown Heights, NY, USA 10598. riolivei@us.ibm.com 
HBIvI Watson Research Center, Yorktown Heights, NY, USA 10598. terhal@watson.ibm.com 
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The work on finding QMA-complete problems was jump-started by a 'quantum Cook-Levin 
Theorem' proved by Kitaev [2] (see also the survey [3]). Kitaev showed that the promise problem 
fc-LOCAL HAMILTONIAN for k = 5 is QMA-complete. Before we state this problem, let us review 
some definitions. A Hamiltonian is a Hermitian operator. A Hamiltonian on n qubits is fc-local 
for constant k if it can be written as Y^j=i Hj where each term Hj acts non-trivially on at most 
k qubits and thus r < poly(n). Furthermore, we require that \\Hj\\ < poly(n) and the entries of 
Hj are specified by poly(n) bits. The smallest eigenvalue of H, sometimes called the 'ground state 
energy' of H, will be denoted as X(H). 

With these definitions in place one can define the promise problem fc-LOCAL HAMILTONIAN as: 

Definition 2 (fc-LOCAL HAMILTONIAN) Given is a k-local Hamiltonian H and a, (3 such that 
(3 — a > po iy( n ) • We have a promise that either X(H) < a or X(H) > (5. The problem is to decide 
whether X(H) < a. When X(H) < a we say we have a 'YES-instance'. 

Kitaev's result was strengthened in Ref. [4], which showed that 3-LOCAL HAMILTONIAN was 
QMA-complete. The subsequent p] proved that also 2-LOCAL HAMILTONIAN is QMA-complete. 

In another direction it was first shown by Aharonov et al. [5] that adiabatic quantum computation 
using 3-local Hamiltonians is computationally equivalent to quantum computation in the circuit 
model. In the adiabatic computation paradigm one starts the computation in the ground-state, 
i.e. the eigenstate with smallest eigenvalue, of some Hamiltonian H(t = 0). The computation 
proceeds by slowly (at a rate at most poly(n)) changing the parameters of the Hamiltonian H(t). 
The adiabatic theorem (see Ref. [6] for an accessible proof thereof) states essentially that if the 
instantaneous Hamiltonian H{t) has a sufficiently large spectral gap, - i.e. the difference between 
the second smallest eigenvalue and the smallest eigenvalue is fi(l/poly(ra))— , then the state at 
time t during the evolution is close to the ground-state of the instantaneous Hamiltonian H(t). 
At the end of the computation (t = T), one measures the qubits in the ground-state of the final 
Hamiltonian H(T). Ref. pQ improved on the result by Aharonov et al. by showing that any 
efficient quantum computation can be efficiently simulated by an adiabatic computation employing 
only 2-local Hamiltonians. 

These results on the complexity of Hamiltonians can be viewed as the first (see also Ref. [7]) in 
a field that is still largely unexplored as compared to the classical case. The class of Hamiltonian 
problems is likely to be a very important class of problems in QMA. Hamiltonians govern the 
dynamics of quantum systems and as such contain all the physically important information about 
a quantum system. The problem of determining properties of the spectrum, in particular the 
ground state (energy) or the low- lying excitations, is a well-known problem for which a variety of 
methods, both numerical and analytical, (see e.g. [HIE]) have been developed. Furthermore, finding 
QMA-complete problems may help us in finding new problems that are in BQP. 

Let us briefly review the classical situation. In some sense the 2-LOCAL HAMILTONIAN problem 
is similar to the MAX-2-SAT problem [10]. But perhaps a better analogue is the set of problems 
defined with 'classical' Hamiltonians such as ISING SPIN GLASS: 

Definition 3 (ISING SPIN GLASS) Given is an interaction graph G = (V, E) with Hamiltonian 

h g = J ij z i TiZi - w 

ijeE i&V 

Here the couplings Jij £ {—1,0,1} and Tj G {—1,0,1} and Z = |0)(0| — |1)(1| is the Pauli Z 
operator. The problem is to decide whether X(Hg) < a for a given a. 

It is known that the problem ISING SPIN GLASS, which is a special case of the 2-local Hamiltonian 
problem, is NP-complete on a planar graph. In fact, it is even NP-complete on a planar graph when 
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Jij = J = 1 and Tj = r = 1 pXJ . In this paper we prove some results on the complexity of a quantum 
version of this model, a quantum spin glass. Our results are based on two ideas. The first one is 
a small modification to the 'quantum Cook-Levin' circuit-to-5-local Hamiltonian construction that 
will prove QMA-completeness of a 5-local Hamiltonian on a 'spatially sparse' hypergraph (to be 
defined below). Such QMA-completeness result on a spatially sparse hypergraph could also have 
been obtained from the 6-dim particle Hamiltonian on a 2D lattice that was constructed in [5]. 

Secondly, we introduce a set of mediator qubit gadget^to manipulate A:-local interactions. These 
gadgets can be used to reduce any fe-local interaction for constant k to a 2-local interaction. Then we 
use the gadgets to reduce a 2-local Hamiltonian on a spatially sparse graph to a 2-local Hamiltonian 
on a planar graph, or alternatively to a 2-local Hamiltonian on a 2D lattice. The general technique 
is based on the idea of perturbation gadgets introduced in Ref. [I]. However the gadgets that we 
introduce here are more general and more powerful than the one in Ref. [Tj. 

Before we state the results, let us give a few more useful definitions. With a 2-local Hamiltonian 
Hq acting on n qubits we can associate an interaction graph G = (V,E) with \V\ = n. For 
every edge in e S E between vertices a and b there is a nonzero 2-local term H e on qubits a and 
b such that H e is not 1-local nor proportional to the identity operator I. We can write Hq = 
X^eS-E He + X^eV where H v is a potential 1-local term on the vertex v. Similarly, with a A;-local 
Hamiltonian one can associate an interaction hypergraph in which the &-local terms correspond to 
hyper-edges in which k vertices are involved. We also use the following definition of a spatially 
sparse hypergraph. A spatially sparse interaction (hyper)graph G is defined as a (hyper)graph in 
which (i) every vertex participates in 0(1) hyper-edges, (ii) there is a straight-line drawing in the 
plane such that every hyper-edge overlaps with 0(1) other hyper-edges and the surface covered by 
every hyper-edge is 0(1). 

A Pauli edge of an interaction graph G is an edge between vertices a and b associated with an 
operator a ab P a (g) P b where P a ,Pb are Pauli matrices X = |0)(1| + |1)(0|, Y = -i|0)(l| + i|l)(0|, 
Z = |0)(0| — |1)(1| and a ab is some real number. For an interaction graph in which every edge is 
a Pauli edge, the degree of a vertex is called its Pauli degree. For such a graph, the X- (resp. Y-, 
resp. Z-) degree of a vertex a is the number of edges with endpoint a for which P a = X (resp. 
P a = Y, resp. P a = Z). 

We will prove the following results. First we show that 

Theorem 4 2-LOCAL HAMILTONIAN on a planar graph with maximum Pauli degree equal to 
3 is QMA- complete. 

With only a little more work, we prove that 

Theorem 5 2-LOCAL HAMILTONIAN with Pauli interactions on a subgraph of the 2-D square 
lattice is QMA- complete. 

Lastly, we answer an open problem in Ref. [5] (see Section [5] for a more detailed statement of 
the result), namely that: 

Theorem 6 Universal quantum computation can be efficiently simulated by a quantum adiabatic 
evolution of qubits interacting on a 2-D square lattice. 

We believe that our Theorem [5] is in some sense the strongest result that one can expect for 
qubits, since we consider it unlikely that 2-LOCAL HAMILTONIAN restricted to a linear chain of 
qubits is QMA-complete. A recent surprising result in this respect is that 2-LOCAL HAMILTONIAN 

*These gadgets are inspired by the idea of superexchange between particles with spin. Loosely speaking, superexchange 
is the creation of an effective spin 'exchange' interaction due to a mediating particle, first calculated by H.A. Kramers 
in 1934 [TJ]. 
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on a one-dimensional lattice with 12-dimensional qudits is QMA-complete |13] . With regards to 
Theorem [U one should note that Aharonov et al. [5] had already proven that interactions of six- 
dimensional particles on a two-dimensional square lattice suffice for universal quantum adiabatic 
computation. Our improvement to qubits on a two-dimensional lattice is an application of our 
perturbation gadgets to [5]'s 6-dim particle construction. 

We would like to draw attention to the power of the perturbative method and in particular to 
the gadgets that we develop in this paper. There are a variety of interesting states that can be 
defined as the ground-states or ground-spaces of k- local Hamiltonians. Prime examples are the 
stabilizer states where the Hamiltonian equals H = I — Y2i Si and S = {Si} is a set of commuting 
stabilizer operators. The ground-space is formed by all states with +1 eigenvalue with respect to 
the stabilizer S and this space is separated by a constant gap from the rest of the spectrum. An 
example is the cluster state [H], the toric code space [15] or any stabilizer code space. Typically, 
the stabilizer operators Si are /c-local with k > 2 which seems to preclude the generation of such 
ground-space as the ground-space of a natural Hamiltonian, see the arguments in Ref. [16]. The 
perturbative gadgets introduced in this paper show how to generate a 2-local Hamiltonian which 
has a ground-space with is approximately a product of a trivial ancilla-qubit space times the ground- 
space of the desired /c-local Hamiltonian. Thus the use of ancillas and the use of approximation 
get us past the constraints derived in [16] . If the original A;-local Hamiltonian has some restricted 
spatial structure, one can show that the resulting 2-local Hamiltonian can be defined on a planar 
graph or, if desired, on a 2-D lattice. 

In the Appendix of this paper we prove a stronger perturbation theorem than what has been 
shown in pQ. The results in the Appendix show that under the appropriate conditions the per- 
turbative method does not only reproduce the eigenvalues of the target Hamiltonian, but also 
the eigenstates, possibly restricted to the low-lying levels of the target Hamiltonian. We believe 
that these results may have applications beyond reductions in QMA and the adiabatic universality 
results in Section [5j 

This paper is organized as follows. In Section [2] we show how to modify Kitaev's original 5-local 
Hamiltonian construction [2j to a 5-local Hamiltonian with interactions restricted to a spatially 
sparse hypergraph. In Section [3] we introduce our perturbation gadgets and in Section 13.11 we 
show how to go from a 5-local to a 2-local Hamiltonian using our basic mediator qubit gadget. In 
Section [3.21 we use new variants of the basic gadget to further reduce the 2-local Hamiltonian on a 
spatially sparse hypergraph to a 2-local Hamiltonian on a planar graph of Pauli degree at most 3, 
Theorem |U With a bit more work we reduce it to a 2-local Hamiltonian on a 2-D square lattice, 
Theorem[5l Finally, Section [5] presents the proof that adiabatic quantum computation using 2-local 
Hamiltonians on a 2D lattice is computationally universal (Theorem [6]) . 

2 A Spatially Sparse 5-local Hamiltonian Problem 

We start by modifying the proof that 5-LOCAL HAMILTONIAN is QMA-complete in Ref. [2] (see 
also [3]). The essential insight is (1) to modify any quantum circuit to one in which any qubit is 
used a constant number of times and (2) make sure that the program to execute the gates in the 
correct time sequence is spatially local. We note that some of the ideas in this section are quite 
similar to those behind the adiabatic 2D-lattice Hamiltonian construction with 6-dim particles in 
Ref. 0. 

Let a quantum circuit use N qubits where n qubits are input qubits and the other N — n qubits 
are ancilla qubits. We first modify this circuit such that gates are executed in R = poly(iV) 
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'rounds' where in every round only 1 (non-trivial) gate is performed |f After a round, the N qubits 
are swapped to a next row of N qubits and then the next gate in the original circuit is executed. 
The total number of qubits in this modified circuit is M = RN. The rows of N qubits for different 
rounds R are depicted in Fig. [U Let us specify an order in which the swap and gate operations 
are executed. In the first round R = 1 we start by applying gates, I and the non-trivial gate, with 
the qubit on the left in Figured! After this round, the swapping starts with the qubit on the right. 
Then again the R = 2 gate-round starts with qubits on the left etc. If we label the gates (including 
/) with a time-index depending on when they are executed, then it is clear that in this model time 
changes in a spatially local fashion. 

We also note that in our construction, each physical qubit enters a gate at most 3 times, twice 
in a swap gate, and once in a I gate or a nontrivial gate. 



R=4 



R=3 



i t i 



R=2 



R=1 



Fig. 1. Two-dimensional spatial layout of the qubits in a quantum circuit for R — 4. A 
qubit is indicated by a •. One and two-qubit gates are indicated by boxes. After the gate 
is executed in row R, those qubits are swapped with the qubits above them in row R + 1. 
The order in which the swap and gate operations are executed can be represented by a 
(time)cursor that snakes over the circuit as follows. We start with the qubit on the left in 
row R = 1. Identity gates are applied on qubits in this row except for the one non-trivial 
gate. We end up at the right and then start swapping the qubits in row 1 with those in 
row 2, starting with the qubit on the right. By doing this we end up at the left. Now we 
perform a round of gate-applications (going right) on the qubits in row R = 2. We end up 
at the right and go left while swapping the qubits in rows R — 2 and R = 3. We continue 
until all necessary gates are executed and the computational qubits are sitting in the last 
row. 



In the class QMA the verifier Arthur uses a verifying quantum circuit V x for an instance x. We 
will use the fact that we can always replace such verifying quantum circuit by a modified verifying 
circuit with the properties that we derived above. 

Given any instance x of a promise problem L 6 QMA and the verification circuit V x , we will 
construct a 5-local Hamiltonian such that 

• if on some input |£,0) V x accepts with probability more than 1 — e (x is a YES-instance), 
then has an eigenvalue less than p ^ n ^ for some polynomial pi(n). 

• if V x accepts with probability less than e then all eigenvalues of are larger than 
for some polynomial P2(n). 

Thus we can map each promise problem in QMA onto a 5-local Hamiltonian problem where the 
specific restricted form of Arthur's verifying circuit will lead to restrictions on the interactions in 



^One could do more gates per round, but this construction is perhaps more easily explained. 
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the 5-local Hamiltonian, that is, the interaction hypergraph will be spatially sparse. In particular, 
when e = 0(2~ n ) for a n qubit proof from Merlin, we obtain a Hamiltonian which obeys the promise 
in Definition 2. Note that Definition 1 uses e = 1/3 but it has been shown, see e.g. [17] , that one 
can make the error e = 0(2~ n ) for a n qubit proof input. 

Thus, these arguments will prove that the 5-local Hamiltonian problem on a so-called spatially 
sparse hypergraph is QMA-hard. Since it is also known that 5-LOCAL HAMILTONIAN is in QMA [2], 
this proves the QMA-completeness of 5-LOCAL HAMILTONIAN on a spatially sparse hypergraph. 

Let us now look at the details of mapping a QMA circuit onto a Hamiltonian problem. Our 
construction is a small modification from the standard construction by Kitaev [2j. We define a 
set of clock-qubits. We use T = (2R — 1)N clock-qubits labeled as c\...,ct- Time t will be 
represented as the state |1*0 T ~*) C1 ... CT as in Ref. [2]. Let U\ . . . Ut be the sequence of operations on 
the computational qubits of the quantum circuit V, one operation for every clock-qubit ci, . . . , ct ■ 
The set of operations includes the actual gates, the / operations when only time advances and the 
swap gates. Let Q; n be the set of n qubits that contain the input Let g ut be the final qubit 
that is measured in the quantum circuit V x . The 5-local Hamiltonian that we associate with 
this circuit is as follows. = H m + -ff ut + -f^ciock + \ J2t=o U cvo \ v (t) where 



H in = |i><i|«®|ioo>(ioo| Ct ,_ 1 , 



Ctq,Ct q + l > 



#out = |0)(0|, out ®|1)(1| CTJ 
T-1 

tfclock = X>l>(0lU lCt+1 . (2) 



i=l 



and 



flevoiv(l) = |00}(00| CliC2 + |10}(10| CliC2 - Lq ® |10)(00| cljC2 - U{ ® |00}(10| C1)C2 , 



^evolv(t) = |100)(100| Ct _ ljCt , C4+1 + |110)(110| Ct _ ljCtiCt+1 

,C t ,C t+ l — U t <8> 1 100) (1 J. v.. I , , ., . ■ 



U t ® |110)(100| ct _ liCtiCt+1 - U\ (8 |100)(110| Ct _ 1)CtiCt+1 , 1< t < T 

t 



^evoiv(T) = |10)(10| CT _ liCT + |ll)(ll| CT _ liCT - U T ® |ll)(10| CT _ ljCT -U T ® |10)(ll| CT _ liCT .(3) 

Hi n is the only term that is different from the 5-local Hamiltonian considered in Ref. [2]; it uses 
the definition of a set of special times t q . Before we define these times, let us look more closely at the 
interactions in the Hamiltonian and how the qubits can be laid out so that each qubit only interacts 
with a set of qubits in its neighborhood. The precise form of this neighborhood is irrelevant, we 
only require that the interaction hyper-graph of this Hamiltonian spatially sparse, as defined in the 
Introduction. 

Given the lay-out of the computational (non-clock) qubits in Figure [U we can 'drape a string' 
of clock qubits over the line following the sequence of computational steps. This ensures that the 
terms in H cvo \ v involve qubits that are in each other's local neighborhood. We can also ensure this 
locality property of H oui by choosing the output qubit g Q ut to be the last qubit on the right in the 
final row. Now let us consider H m . For every qubit in the layout in Figure [T] there is a time in 
which the running cursor which snakes over the circuit first arrives at this qubit. For the qubits 
in R = 1, this is when the cursor comes from the left doing the I operations or the non-trivial 
gate. For the qubits in the other rows R > 1, it is when the cursor, coming from the right, starts 
swapping the qubit with the previous row R — 1. These cursor actions are represented in H em \ v . 
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For a qubit q we define the clock-qubit c% as the clock-qubit whose bit is flipped in the interaction 
representing the earliest gate (the action of the cursor) on the qubit q in H evo \ v . Then it is clear 
that the clock-qubit ct q is local to the qubit q and therefore H 1Ti again represents an interaction 
between qubits that are in each other's local neighborhood. It is also clear that the role of H m is 
to make sure that the state of the qubits is set to before the gates actually act on these qubits. 
Note that we set the state of all qubits (except those in Qin) to zero, also the ones in the later rows 
that are merely used as dummy qubits to be used in swaps. This is not absolutely necessary but 
merely convenient. 

These arguments show that the interaction hypergraph of the Hamiltonian is spatially sparse. 
Note also that given a quantum circuit with N qubits one can efficiently construct the interac- 
tion hypergraph of the corresponding Hamiltonian and draw this hypergraph in the plane where 
hyperedges involving 5 qubits are represented as five-sided polygons. 

The proof of the following Lemma is analogous to the proof of Theorem 14.3 in [2]. 

Lemma 1 Let \f) = yf^Y^a l&) ?1 ... 3 Jl t T -*) Cl ... CT where |&) = U t \&-l) for all 1 < t < T 
and |£o) = n ) for some state |£) of the input qubits. If Arthur's verifying quantum circuit 

V x accepts with probability more than 1 — e on some input 00 ... 0) then (^\H^\^) < IfV x 
accepts with probability less than e on all inputs |£,0) then all eigenvalues of are larger than 
or equal to c<yl ^ — for some constant c. 

Proof. Consider first (ip\H^ We only need to check that {^\H- m \^) = since this term is 
different than the one in Ref. [2]. We note that H in \^} oc E^Q in I 1 )( 1 lgl6 q -i, l t »- 1 T- (*«- 1 )) = 
since in all computational qubits are set to before they are being acted upon, i.e. qubit q is the 
state at all times t < t q . Thus \if)) has zero eigenvalue with respect to all terms in except H out . 
If V x accepts with probability more than 1 — e, this implies that (ip\H^\ip) = (ip\H on t\ip) < jqr[- 
The second part of the proof is to show that if V x accepts with small probability, the eigenvalues 
of H are bounded from below. Again the proof is identical in structure to the proof in [2] except 
for H m . We first note that preserves the space of 'legal' clock-states 5, i.e. clock-states of 
the form |l*0 T_t ) and thus we can consider the minimum eigenvalue problem of on S and S 1 - 
separately. On S 1 - this minimum eigenvalue is 1 since at least one of the constraints of -£f c iock is 
not satisfied. Now we consider H^\s which we can express using the definition \t) = |1*0 T ~*). We 
have Hi n \s = Yq£Q in \^){Mq ® — !)(*<? — 1|- As in the standard proof we perform a rotation W 
to a more convenient basis where W = Ylt=o Ut ■ ■ - U\® \t){t\. Let 

H 2 = H cvolv \ s W = I®E, (4) 

where E is defined below Eq. (14.9) in [2]. Let 

H x = W\H m + H out )\ s W = |1><1|« ® \h ~ !)<*? - 11 + C/ t |0)(0| 9out C/ ® \T)(T\, (5) 

q^Qin 

where U = Ut ■ ■ - U\. Note that H- m \$ is unchanged by the rotation W since there are no gates 
acting on a qubit q prior to the time t q . Now we would like to use Lemma 14.4 in Ref. [2] and 
lower-bound the smallest eigenvalue of H1 + H2. Let C\ and £2 be the non-empty null-spaces of H\ 
and H2- Lemma 14.4 states that for such H\ > and H2 > Owe can bound H 1 + H 2 > 2v sin 2 (6/2) 
where v is the smallest non-zero eigenvalue of H% and H2 and cos 2 6 = max, 6 £ 2 (rj\Pci \v) where 
P£j is the projector on C\. The minimum of the smallest non-zero eigenvalue of H\ and H2 is as 
in Ref. [2], namely v > cT~ 2 . 
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Now we show that, as in [2], one can bound sin 2 6 > 1 ■ Putting these results together 
shows that the minimum eigenvalue of is at least c<yl ~^^ for some constant c, as claimed. 
As in Ref. [2] any state in £2 is of the form |£) ® "^fqrf Ylt=o I*) where |£) is arbitrary. We can 

also write P A = ^=0 p t ® |*> <*| where P T = rf\l){l\ qout U, and P t = n^ Qin | t(J=m |0)(0| ? tg> / cW 
where / e lset is the / operator on all computational qubits for which t q 7^ t + 1. At some times Pt 
may just be / on all qubits. Here R q dQ in \t =t+i is tensor product of |0)(0| for all qubits q for which 
tq = t + 1. Thus we need to bound 

cos 2 # = ^max(£|^Pi|0. (6) 

All Pt for t < T commute and their common eigenspace is the space where all qubits q ^ Q m are 
set to 1 00 . . . 0). We can write any |£) as |£) = a|00 ... 0, ?/>o) + \f3) where ipo is a state for all qubits 
in Q m and \(5) is a state with norm 1 — |q| 2 in which at least one of the k non- input qubits is not 
in |0). Thus we have 

cos 2 < (|a| 2 T + |a| 2 (0^o|PT|0,^o) + 2|a| \(0,^ \P T \P)\ + (T-l)(0\/3) + (P\P T \P)) • (7) 

Given the acceptance probability of the circuit V x we can bound (0, iPo\Pt\0, tpo) < e. We also 
bound (f3\P T \f3) < (I3\0). This gives 

cos 2 9 < ^ (T + \a\ 2 e + 2\a\ Jl^W) < 1 - 1 ~ "'f . (8) 

□ . 



3 Perturbation Theory 

In this section we introduce the perturbation method. Our main new idea is the use of mediator 
qubits that perturbatively generate interactions. The mediator qubits are weakly coupled to the 
other qubits and to lowest order in the perturbation this coupling generates an interaction between 
the other qubits, see Section 13. 11 We will show as a first step how this can be used to reduce any 
/c-local Hamiltonian problem to a 3-local Hamiltonian problem. We can then use the perturbation 
gadget in pQ to reduce a 3-local to a 2-local Hamiltonian (we also sketch an alternative mediator 
qubit method) . To reduce a 2-local Hamiltonian to a 2-local Hamiltonian on a 2D lattice or a planar 
graph, we need a few other applications of our mediator qubit gadgets which will be introduced in 
Section 13.21 

In Ref. pQ the authors reduce the problem 3-LOCAL HAMILTONIAN to 2-LOCAL HAMILTONIAN 
by introducing a perturbation gadget. The idea is to approximate A(P targ et) of a desired (3-local) 
Hamiltonian -fftarget by \(H) of a 2-local Hamiltonian H where \(H) is calculated using pertur- 
bation theory. One sets H = H + V where H is the 'unperturbed' Hamiltonian which has a large 
spectral gap A and V is a small perturbation operator. We will choose H such that it has a de- 
generate ground-space associated with eigenvalue and the eigenvalues of the 'excited' eigenstates 
are at least A. The effect of the perturbation V is to lift the degeneracy in the ground-space and 
create the target Hamiltonian in this space. 

More accurately, we have a Hilbert space C = where £_ is the ground-space of H. Let n.|- 

be the projectors on C±. For some operator X we define X ++ = n + Xn + ,X |_ = n_Jn + ,X-| = 

n_|_ATI_, X = n_ATI_ and X + = X ++ . In order to calculate the perturbed eigenvalues, one 

introduces the self-energy operator S_(z) for real- valued z 

E-(z) = H- + V— + V- + G + (I + - V ++ G + )- 1 V+-, (9) 
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where we can perturbatively expand 

(/+ - V++G+)- 1 =I+ + V ++ G + + V ++ G + V ++ G + + .... (10) 

Here G+, called the unperturbed Green's function (or resolvent) in the physics literature, is defined 
by 

G+ 1 = zl+ - H+. (11) 

In Ref. [T] the following theorem is proved (here we state the case where the ground-space of H 
has eigenvalue and H has a spectral gap A above the ground-space): 

Theorem 7 ([1]) Let \\V\\ < A/2 where A is the spectral gap of H and X{H) = 0. Let H\ <l \/2 
be the restriction of H = H + V to the space of eigenstates with eigenvalues less than A/2. Let 
there be an effective Hamiltonian H e g with Spec(i/ C fj) C [a, b]. If the self-energy H-(z) for all 
z £ [a — e, b + e] where a < b < A/2 — e for some e > 0, has the property that 

\\E-(z) - Hrf\ \ < €, (12) 

then each eigenvalue Xj of H\ < /±j 2 is e-close to the jth eigenvalue of -ff e ff- In particular 

\X(H eS ) - X(H)\ <e. (13) 



This theorem can be generalized to Theorem I A. II proved in the Appendix. Theorem I A . 1 1 shows 
that under appropriate conditions, the effective Hamiltonian is approximately identical to H re- 
stricted to its low-lying eigenspaces. With the same technique we also prove Lemma I A. II in the 
Appendix which shows that the ground-space of a target Hamiltonian can be generated perturba- 
tively (under the assumption that the target Hamiltonian has a l/poly(n) gap). Lemma I A. II was 
also proved in [1] in the special case that the ground-space is non-degenerate. 

3.1 Mediator Qubit Gadgets 

In the following explanation of the gadgets we will refer to -fftarget as the desired Hamiltonian that 
we want to generate perturbatively and the effective Hamiltonian is H e R = iif target ® |00 . . .) (00 . . . | , 
i.e. the ancillary 'mediator' qubits are in their ground-state 1 00 ... 0). 

The gadgets that we introduce below to accomplish the reduction are what we call mediator 
qubit gadgets and seem to be useful in general to manipulate fc-local interactions. The idea is 
that we replace a direct interaction between two groups of \k/2] qubits with indirect interactions 
through a mediator qubit. In the ground-state of the unperturbed Hamiltonian H the mediator 
qubit is in state |0). The perturbation V is chosen such that interaction with the other qubits 
can flip the mediator qubit. The perturbative corrections to the self-energy, up to second order in 
the perturbation, involve the process of flipping the mediator qubit by interaction with a group of 
qubits a and flipping the mediator qubit back to |0) by a second interaction with a group of qubits 
b. If a = b we potentially obtain some \k/2] -local terms. For a^&we obtain an effective fc-local 
interaction involving groups a and b. This gadget could also be used with three or more groups 
of qubits (or higher dimensional quantum systems); in this case interactions would be generated 
between all groups of qubits. An example of such application is the Cross gadget, explained in 
Section ET21 

Subdivision Gadget. Assume that a fc-local operator associated with (hyper)edge ab is of the 
form A (g> B and let r = max(||^4||, ||5||). The hyper-edge ab is part of a larger (hyper)graph and 
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A B A X X B 
• • • o • 

a b a w b 

Fig. 2. Subdivision gadget. A fc-local interaction is reduced to \k/2] + 1-local interactions 
using a mediator qubit vertex w. The operators A,B,X next to the edges indicate which 
operators correspond to the edges. 



a corresponding Hamiltonian. Let all other terms in the Hamiltonian be H e \ se . We can write the 
Hamiltonian as 

^target = (#cise + A 2 /2 + B 2 /2) - {-A + B) 2 /2 = fl^, - (-A + £) 2 /2, (14) 

so that -ffg lse contains some additional \k/2] -local terms as compared to H e \ sc . W.l.o.g. we assume 
that max(||i^ lse ||,r) > 1. 

The terms in the gadget Hamiltonian H = H + V are the following 
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ff = A|l)(l| ro , V = H' else + y/A/2(-A + B)®X w . (15) 

The operator X w is the Pauli X operator acting on qubit w. The degenerate ground-space £_ 
of H has the mediator qubit in the state |0). We have the following: iT_ = 0, G+(z) — 
V—=H f e]se ®\0)(0\ u] and 

V+- = y/A/2(-A + B)®\l)(0\ vl . (16) 

Thus the self-energy equals 

E-W = (^cisc + ("A + i?) 2 ) |0)(0U + O ((j^a) ■ ( 17 ) 

We can expand the self-energy around z = and identify H c q = ^target ® |0)(0|. This gives 

l|s _ w _ ffeffll . (^!) + o(W) +0 (NW). (18) 

In order for Theorem [7] to apply the following must hold: (1) for z G [— ||-ff e ff|| — ||i? e fr|| + e ]> 
should be e-close to H e g and (2) ||V|| < A/2. Let us consider how to choose A such that 
these conditions are fulfilled. We can bound ||V|| < ||#e lse || + ^2A r ^ (\\ H ' e \se\\ + V2r). We 
will choose A such that \z\ < A. Then, using the bound on ||V|| gives 

r 2 \ „ { (\\H' cl J\ + V2rf\ 



- H eS \\ <0^-j+0 ^ ' j . (19) 

Let us choose 

A=(||^ lsc ||+C 2 r) 6 /6 2 , (20) 

for some constant C2 > V%- This choice lets us bound the last term in Eq. (fl9|) by 0(e). Since 
A- 1 < A , we can bound the first term in Eq. (|19p by 0(e 2 )- Let us verify the second condition 
< A/2 with this choice of A. We have indeed 

< = — < e. (21) 

" (11^,11 + V2r) 2 " 
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Consider the conditions on \z\, i.e. z £ [— ||Heff|| — e, ||i? e ff|| + e] and \z\ < A. Since ||i? e fr|| < 
1 1 i^ci se 1 1 + 2r 2 , we can consider the interval \z\ < \\H' eise \ \ + 2r 2 + e. For sufficiently small e we have 
(using max(||ffg lse ||,r) > 1) 



\z\ e 



2/ 



^clsell+2r 2 +6) 



A (||^|| + C 2 r 



< 0(e 2 ) < 1. (22) 



Thus for the choice of A as in Eq. (fl9j) we have S_(z) = ^target ® |0)(0| w + O(e). From Theorem 
[7] it follows that \X(H e s) — X(H)\ = O(e). When ||i7g lse ||, r and 1/e are polynomial in n (n is the 
number of qubits of i?target)> h is clear from Eq. (|20p . that the norm of the gadget Hamiltonian 
H which uses A is polynomially larger than the norm of the effective Hamiltonian. This implies 
that the gadget can only be used a constant number of times in series if norms have to remain 
polynomial. 

We will use this type of gadget in parallel, that is, in many places in an interaction graph 
at once. Let us explain how this happens in detail and argue that the local gadgets operate 
independently, i.e. there are no cross-gadget contributions to 2nd order in the perturbation. Let 
^target = #eise - Yn=i ^target where #t ar g e t = (-A* + Bi) 2 /2 for some operators A* and Bj. H else 
contains all interactions that are not generated perturbatively in addition to the compensating terms 
Af/2 etc., similar as above. We introduce k mediator qubits W\ . . . Wk and choose H = J2i Hi + V 
where iJ; = A|l)(l| Wi and V = H clsc + ^Afi^-Ai + B<) ® X Wi . 

The degenerate ground-space of H has all mediator qubits w\ . . . in the state |0). Let h(x) 
be the Hamming weight of a bit-string x G {0, l} k of the qubits w\ . . . w^. We have the following: 
G+ = E^oo...o T^K^ V— = H else ® |00. . . 0)(00 . . . 0| and 

y+_ = yAT2^(-A i + S i )|00...1 i ...0)(00...0|, (23) 

i 

where 1 00 ... lj ... 0) has qubit Wi in the state |1). To second order in the perturbation V, there 
are no cross-gadget terms in S_(z). Thus the self-energy to second order equals 

S_ (z) = (tf clsc + 2(z A A) ^ {-Ai + Bif^ ® |00. . .0)(00. . .0| +0 ( J^aV ) ' (24) 

Choosing A = poly(n)/e 2 for some sufficiently large poly(n) gives 

E_(Z) = target ® |00 . . . 0) (00 . . . 0| + O(e). (25) 

We need to use the parallel application of this gadget twice in order to reduce the ground- 
state energy problem of our 5-local Hamiltonian to that of a 3-local Hamiltonian; one application 
results in a 4-local Hamiltonian, another one reduces it to 3. Similarly, any fc-local Hamiltonian for 
constant k can be reduced to a 3-local Hamiltonian by these means. A 3-to-2-local reduction can 
be carried out using the gadget in pp. However an alternative construction exists which we now 
explain. 

3-TO-2-LOCAL gadget. Assume that we have a target Hamiltonian //target = A B C + i/eise- 
The idea is to generate the 3-local term A ®B®C by using perturbative effects up to third order. As 
before one introduces a mediator qubit w whose ground-state is |0) for the unperturbed operator. 
And, as before, we have perturbations proportional to A <g> X w and B ® X w which can flip the 
mediator qubit. We also have a perturbation V which contains a term proportional to C <S> 
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which implies that there is an interaction with C if the mediator qubit is 'excited'. Thus, the 
second-order perturbative corrections give us terms proportional to A <g> B whereas third-order 
corrections gives us the desired A® B <g> C (and some additional 2-local terms). More precisely, let 
^target = -^eise + A (g> B <g> C. Let r = max(| \A\ I , I \B\ 1 , 1 1 C\ I) . We choose H = A|1}(1|«, and 

V = H e]Be + F extra - A 2 / 3 C ® + A 2 / 3 (-A + B) ® X w /V2 (26) 

where the additional 2-local compensating term is V ex tra = A 1 / 3 (—A + B) 2 /2 + (A 2 + B 2 ) ® C/2. 
One can show that 

£_(*) = [F else + A ® B ® C] ® |0)(0U + O(|z|A- 2 / 3 ) + O(A- 1 / 3 ). (27) 

For sufficiently large A and |z| < ||i? e ise|| + 0(r 3 ) + e we make 'E-(z) sufficiently close to -f/target <8> 
|0)(0|. 

The important conclusion of this section is that one can derive a 2-local Hamiltonian on a 
spatially sparse graph for which the ground-state energy problem is QMA-complete. The interaction 
graph is restricted because the perturbation gadgets preserve the spatial restrictions of the original 
hypergraph of the 5-local Hamiltonian. 

3.2 More Mediator Qubit Gadgetry 

For our next round of reductions we need to describe some different uses of the subdivision gadget 
acting on 2-local interactions. In the following we will assume that every edge in the interaction 
graph is a Pauli edge. It may thus be that the interaction graph contains other edges between the 
same vertices, each edge associated with a different product of Paulis. The Pauli degree of a vertex 
is then the number of Pauli edges that are incident on this vertex. 

The Cross Gadget. For the Cross Gadget we assume that we have a graph G which, when 
embedded in the plane, contains two crossing edges such as in Fig. El Assume that the operator 
on edge ad is a ad P a ® Pd an d on edge be we have a bc P b <8> P c - Our desired Hamiltonian is 

^target = #else " ("CW^a " a bc P b + P c + P d f /2. (28) 

It is clear that the last term in this Hamiltonian generates the desired crossing edges a a( iP a <g> P^ 
and ab c Pb <S> P c hi addition to other operators on the edges ab, bd, cd and ac. Thus -ff e i se is a sum of 
all other operators associated with the original graph G and a set of operators on the edges around 
the cross, see Figure El which are meant to cancel the extra operators generated by the last term 
in Htarget- As before we set H = H + V with 

# = A|1)(1U V = H else + y r A/2(-a ad P a -a bc P b + P c + P d )®X w , (29) 

and the analysis follows as for the subdivision gadget. Note that if there are no edges ab, bd, cd, 
or ac in iJtarget; there will be such edges in H, as indicated in Fig. [3j 

The Fork Gadget. For the Fork gadget we have a subgraph as in Fig. |4] where the operator 
on edge ab is a ab P a ®P b and on edge ac it is a ac P a ®P c . The Fork gadget merges the 2 edges coming 
from vertex a at the cost of creating an additional edge between b and c. Our desired Hamiltonian 
is 

^target = #else ~ {Pa ~ Oi ab P b - a ac P c ) 2 /2, (30) 

where -ff e ise contains all other terms not involving edge ab and ac. We take 

H = A\1)(1\ W , V = H elBe + ^/A/2(P a -a ab P b -a ac P c )®X w , (31) 
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Pd 



Fig. 3. Cross gadget. A crossing between two edges is removed by placing a mediator 
qubit in the middle. Additional edges ab, ac, bd and cd are created. 




Fig. 4. Fork gadget. Two edges of the same type at vertex a are merged by the placement 
of a mediator qubit w. The additional edge be is created. 

and the analysis follows as before. 

The Triangle Gadget The Fork gadget can also be used in order to reduce the degree of 
a vertex, see Fig. [SJ this is achieved by applying the Fork gadget together with the subdivision 
gadget in series. We first apply a subdivision gadget on the edges ab and ac. Then we apply the 
Fork gadget on vertex a, thus generating the inner triangle in Fig. 




Fig. 5. Triangle Gadget. We first subdivide edges ab and ac and then apply the Fork 
gadget on vertex a. This give rise to a 'mediator triangle' such that vertices b and c have 
the same degree as before and vertex a has reduced its degree by 1. 

4 2-LOCAL HAMILTONIAN on a 2-D Square Lattice 

With these tools in place, we are ready to state the reduction which we obtain by applying the 
gadgets in the previous section. Together with our previously argued 5-local to 2-local reduction, 
this Lemma implies Theorem HI 

Lemma 2 Let Hq be a 2-local Hamiltonian related to a spatially sparse graph G = (V, E) with 
\V\ = n and where Hq = ^2 e&E H e + YlveV H v such that \\H e \\ < poly(n) and \\H V \\ < poly(n). 



13 



Fig. 6. Localizing a vertex. 

For any e > there exists a graph C sim which is planar with maximum Pauli degree at most 3 and 
a polynomially bounded 2-local Hamiltonian H Gs im such that 

\\(H G )-X(H Gsim )\=0(e). (32) 

Moreover, there is a planar straight-line drawing of G sim such that all edges in G sim have length 
0(1), and all angles between adjacent edges are 0(1). 

Proof. 

• We use the subdivision gadget in order to localize each vertex with Pauli degree more than 
3, see Fig. [6j Then we are ready to reduce the Pauli degree (which is some constant) of these 
vertices. 

• Consider the set of vertices with Pauli degree more than 3. We are going to apply the Triangle 
gadget to all these high degree vertices in the following way. We first apply the subdivision 
gadget to all edges that we intend to merge using the Fork gadget; we can do this in one 
parallel application. We do this so that the triangle gadgets that we will apply in parallel 
never act on the same edges. Then, for a vertex with X-degree d x , Y-degree d y , Z-degree 
d z we do the following. We pair the X-edges and apply to each pairing a Fork gadget. This 
means we have reduced the X-degree to \d x /2\. In parallel we pair the Y-edges and the 
Z-edges using the Fork gadget, halving their degrees. We do this single perturbative step in 
parallel for all high-degree vertices in the graph. We repeat this Triangle gadget process 0(1) 
number of times (since the maximum degree initially was O(l)) until the total Pauli degree 
of every vertex is at most 3. Since the initial degree of every vertex was O(l), the number of 
additional crossings that we generate per edge is constant. 

• Next, we reduce the number of crossings per edge, by subdividing each edge a constant number 
of times, see Fig. UJ Every subdivision is done in parallel on all edges of the graph that need 
subdividing. 

• Then we use the subdivision gadget to localize each crossing, see Fig. We apply the 
subdivision gadget in parallel on every crossing in the graph and we repeat the process 4 
times so that for all crossing edges ab, cd, the quadrilateral acbd contains only these points 
and the crossing edges. 

• We apply the Cross gadget, see Fig. O in parallel to every localized crossing in order to 
remove the crossing. Note that due to the localization step the cross-gadget only involves 
mediator qubit vertices with degree at most 2. Thus the cross- gadget generates additional 
2-local terms around the square, but the total Pauli degree of the resulting vertices is at most 
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4. Note that these vertices with degree 4 are all mediator qubits which only have non-zero 
X-degree (and zero Y- and Z-degree). 

• On all mediator qubits with X-degree 4 we apply the Triangle gadget reducing the degree to 
3. Since the triangle gadget generates mediator qubits with X-degree 3 we cannot do any 
further reductions. 

Thus in this final Hamiltonian there are no vertices with Pauli degree more than 3 and the graph 
is planar. Theorem \7\ is used in every gadget application to give the final result, Eq. (|32p . Note 
that by this reduction all original system qubits have Pauli degree at most 3 by having X-degree, 
Y-degree and Z-degree ranging from to 1. The mediator qubits have X-degree ranging from 2 to 
3 and Y- and Z-degree. □. 



• • • o — o — o — o • 

Fig. 7. An edge that crosses C other edges is subdivided [log C] times by inserting a 
mediator qubit. 



b 




Fig. 8. Localizing a crossing by applying the subdivision gadget four times. 



4-1 Representation on a 2-D Square Lattice 

Any planar graph G = (V, E) with maximal degree 3 in which the (straight-line) edges have length 
O(l) and adjacent edges form an angle of fi(l) can be represented on a planar square lattice in 
the following sense: each vertex a of G is mapped to some lattice site 4>(a) inside the square 
[— 0(|V|), 0(|F|)] 2 , and each edge ab of G is mapped to a lattice path (f>(ab) of length 0(1) from 
4>{a) to (f)(b) that does not cross any other vertices or any other path. To see this, one can look at 
Fig. [U]or follow these steps: draw a fine square grid on the plane. If the spacing between points on 
the grid is small enough, moving each vertex a of G to a vertex in the lattice (and redrawing the 
edges) still leaves the graph planar, with 0(l)-length edges and £1(1) angles. Now for each edge, 
draw a lattice path that stays close to the edge. If the grid is fine enough, these paths can never 
cross outside an 0(l)-size square (indicated in grey in Fig. [9|) around the vertices of the graph, 
because of the angle condition. By further refining the grid if necessary, one can reroute each of 
the paths stemming out of a vertex a inside of a's square, so that no two different paths collide. It 
is easy to see that we only need the grid to have spacing 0,(1), and that all the other conditions 
above are satisfied. 

Clearly, this embedding can be found efficiently, given the adequate embedding of G. If H is a 
Hamiltonian that has G as (Pauli) interaction graph, one can use the subdivision gadget 0(1) times 
in parallel to map each edge ab to a path of the same length as <j)(ab). The Hamiltonian H thus 
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Fig. 9. A planar graph of maximal degree < 3 and its representation in the lattice. In the 
gray squares, the paths are rerouted to avoid crossings. 

obtained has interaction graph 4>{G) and X(H) is 0(e)-close to X(H). These arguments together 
with our previous results and Lemma [2] prove Theorem [5j 

5 Universal Quantum Adiabatic Computation 

In Ref. pQ the authors show that their perturbation-theoretic reduction of 3-LOCAL HAMILTONIAN 
to 2-LOCAL HAMILTONIAN also reduces 3-local adiabatic computations to 2-local ones. The goal 
of this Section is to show that an analogous result can be carried out in the present context, namely 
that 2-local Hamiltonians with nearest-neighbor interactions on qubits on a 2D lattice suffice for 
universal adiabatic quantum computation. 

Let us describe in more detail what our goal is. We will construct a (classically) poly-time 
computable map $ that takes as input a classical description (Q) of a quantum circuit Q and 
outputs a description of an adiabatic quantum computation on a 2D lattice. Suppose Q acts on n 
qubits and has T gates. Then 

Here 

1. p £ N is a constant independent of Q; 

2. N, the number of qubits on which Hi acts is poly(n,T); 

3. for each i € {0, . . . ,p}, (Hi) describes a 2-local nearest-neighbor Hamiltonian on qubits, acting 
on the same subset of N = poly(n, T) sites of the square lattice; 

4. \\Hi\\ = poly(n,T) for all i G {0, . . . ,p}; 
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5. Let 

v 

^ = ^^,(^[0,1]). 

i=0 

The spectral gap between the ground-state and first excited state of H(s) is l/poly(n, T) for 
all s. 

6. The ground-state of H(0) is |0)® and the ground-state of H(l) encodes the result of the 
computation of Q on input |0)® n (a more precise description is given in [H Section 7] or [5]). 

Of course, all occurrences of poly above correspond to fixed polynomials that do not depend on 
Q. Notice that for any Hamiltonian satisfying the above conditions one has that 

< poly(n,T), for all j = 0, 1, . . . 

This is sufficient to ensure that adiabatic computation implemented by H(s), starting from |0)® , 
appropriately simulates the quantum circuit Q, that is, in polynomial time [6j. Note that the usual 
adiabatic computation, e.g. the universal adiabatic computation in [5], has p = 1. 

There are several ways to map a circuit Q to a corresponding H(s). One could modify the 5-local 
construction in this paper in order to show that one can do universal quantum computation using 
a quantum adiabatic computation with a 5-local Hamiltonian on a spatially sparse graph. Then 
we could apply the perturbation gadgets to derive a 2-local Hamiltonian with similar properties. 
However, an easier route to the desired result is the following. In [5j it was shown how to map a 
circuit Q to a corresponding Hamiltonian H^(s) on a 2D lattice. That construction satisfies all 
but one of the above requirements, as it acts on 6-dimensional qudits rather than qubits. However, 
we can embed 6-dimensional qudits in states of 3 qubits. This implies that the 2-local interac- 
tions between these particles will be mapped onto 6-local interactions. Then we can apply the 
perturbation gadgets to 'massage' this Hamiltonian on a spatially sparse hypergraph to a 2-local 
Hamiltonian as we have done in our QMA construction. 



sup 

sefo.il 



6PH(a) 



Let us first review the 6-dim particle Hamiltonian, see Sec. 4.2 in [5]. The four phases of the 
particles, the unborn, the first, second and dead phase, can be described by two qubits in the 
states |unborn) = |00),|first) = 1 01) , (second) = 1 10) and |dead) = |11). The third qubit holds 
the actual computational degree of freedom. In the 6-dimensional representation of the unborn 
and dead phase the computational degree of freedom is assumed to be fixed. If we represent those 
states as 3-qubit states, we fix the third qubit to be in the state |0). Hence we obtain 6 states: 2 
'first' states 1 010) , 1 011) , two 'second' states 1 100) , 1 101) and one unborn state 1 000) and one dead 
state 1 110). With this mapping the entire Hamiltonian in Sec. 4.2 in [5] can be rewritten in terms 
of 6-local interaction between qubits. Since we embed the 6-dim particle Hamiltonian in a higher 
dimensional space, we need to make sure that states outside the embedded space (i.e. |001) and 
| 111)) are penalized in the Hamiltonian, i.e. do not contribute to the ground-space. In Table 1 in 
[5] a list of forbidden configurations is given. In this list we can replace every unborn state by two 
unborn states | unborn, 1) and |unborn, 0) and similarly for the dead states. This implies a small 
modification of H" ioc ^- As a consequence we get that the space of legal shapes S is the same for 
this embedded Hamiltonian as for the original 6-dim particle Hamiltonian. It then follows that one 
can apply Lemma 4.6 and 4.7 bounding the spectral gap of the Hamiltonian in the space of legal 
states. We note in passing the Hamiltonian (s) has only linear terms in s (that correspond to 
p = 1 above) and terms independent of s (corresponding to p = 0). 
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Our second step is to analyze how this desired 6-local Hamiltonian can be implemented using 
a 2-local Hamiltonian on a 2D lattice. It is clear that one can apply the perturbation gadgets in 
Section 3 and 4 of this paper and map a 6-local Hamiltonian on a spatially sparse hypergraph onto a 
2-local Hamiltonian on a subgraph of the 2D lattice. To go from a 3-local to a 2-local Hamiltonian 
we will use our alternative 3-to- 2-local gadget described in Section 13.11 One needs to show the 
following properties of the perturbation method in order for these reductions to work: 

1. The 2-local adiabatic path Hamiltonian H^ 2 \s) obtained through the perturbation gadgets 
simulates the 6-local adiabatic path Hamiltonian. This implies that the ground-state of the 
2-local Hamiltonian should be approximately the ground-state of the desired 6-local Hamil- 
tonian and the gap for the 2-local Hamiltonian is approximately the gap of the 6-local Hamil- 
tonian. This requires showing that the perturbative method that we employ does not only 
reproduce the lowest-eigenvalue but also the ground-state and the gap above the ground-state. 

2. One needs to verify that H^ 2 \s) is of the form Yli=o s * Hi, with p constant and maxj \\Hi\\ < 
poly(n,T). 

Our 2-local simulator Hamiltonian H^ 2 \s) is determined by applying the perturbative gadgets 
in Sections 13.11 and 13.21 on the 6-local Hamiltonian H^(s). In [1] it was shown how to generate, 
not only the lowest eigenvalues, but also the ground-state with the perturbative technique. This 
implies that both the ground-state of the target Hamiltonian as well as the gap above this ground- 
state can be generated perturbatively. Since the total number of applications of the perturbation 
theory is constant, one can apply this argument for each step and thus show that the 6-local target 
Hamiltonian can be effectively generated by a simulator Hamiltonian . 

We now fulfill our second task, i.e. we show that H^ 2 \s) = Yl\=o s% ^i, with p constant and 
\\Hi\\ polynomially bounded. In [1] such arguments were developed for the 3-to-2 local perturbation 
gadget and basically identical arguments can be given here. The original Hamiltonian is at 
most linear in s. If a gadget is applied on a term which is linear in s, for example a 6-local term 
such sA®B = A(s) ® B, we obtain a new Hamiltonian of which the terms are at most quadratic in 
s. Similarly each application of the perturbation gadgets takes a Hamiltonian H'(s) = Yl^=o sJ Hj 

to another Hamiltonian H"(s) = Yli=o s * H'( where p" < 2p' . Assuming that the norm of each H[ is 
polynomial in n and T, then the norms of each H'J are also poly(n,T). Thus the final Hamiltonian 
H( 2 \ obtained after a constant number of gadget applications, is indeed of the desired form. 

6 Discussion and Acknowledgements 

The drawback of the reductions performed by our perturbation theory method is that the 2-local 
Hamiltonian that we construct has large variability in the norms of the 2-local terms. In other 
words, 2-local terms have constant norm whereas others can be fairly high degree polynomials in 
n. Such dependence on n may be undesirable from a practical point of view, e.g. if one wants to 
perform universal adiabatic quantum computation. 

It is possible that a less stringent but still rigorous perturbation theory could be developed 
in which only the expectation values of local observables with respect to the ground-space are 
perturbatively generated. If such expectation values are reproduced with constant accuracy (not 
scaling as l/poly(n)), then the perturbation theory need not be accurately reproduce the entire 
ground-space as in Lemma 3. For adiabatic quantum computation this method would suffice since 
one can measure a single output qubit to extract the answer of the computation. 

One of the reasons why finding QMA-complete problems is of interest is that it may give us 
a hint at what problems can be solved in BQP. One example is the unresolved status of the 
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2-local Hamiltonian problem on qubits in one dimension. Another example is to find a quantum 
extension of classical 2-local Hamiltonian problems which can be solved efficiently. We thank David 
DiVincenzo for an inspiring discussion about super exchange. We would like to thank Sergey Bravyi 
for pointing out an improvement in the proof of Lemma 2. We acknowledge support by the NSA 
and the ARDA through ARO contract number W911NF-04-C-0098. 



Appendix A 



7 General Perturbation Theorem 

In order to give a more complete background in the perturbation method we will prove in Theorem 
IA.1I that under the right conditions the entire operator i?|<A» is approximated by H e f[, not only 
its eigenvalues. In Ref. [1] a similar result was proven, namely that the ground-state of H is 
approximately the ground-state of H CT[ . We extend their result to the case when the ground-space 
is degenerate in Lemma I A. II of this Appendix. To a certain extent our proof-technique is similar 
to the one used in Ref. pQ , however we will use complex z and contour integration in parts of the 
proofs. 

Some of our notation has been given in Section [3] for the specific cases considered in this paper. 
Here we consider the more general setting as defined in Ref. pp. 

Assume that H and V are operators acting on the Hilbert space C and H = H + V . H has 
a spectral gap A such that no eigenvalues lie in the interval [A* — A/2, A* + A/2] for some cutoff 
A*. Let £_ (resp. C + ) be the span of all eigenvectors of H whose eigenvalues are less than A* 
(respectively larger than A*). We will use the resolvent G(z) = (zl — H)^ 1 of H with complex 
z £ C and let G(z) = (zl — H)~ l be the resolvent of H. The definition of the self-energy is 
given by 

£_(*) = zl- - GZl(z). (A.l) 



see also Eqs. (|9|)-(|T0|). 

The perturbation theory result of Kempe et al. states that under suitable technical conditions, 
-namely if X_(z) is close to a fixed operator H e $ for all z in some range-, all eigenvalues of 
H = H + V that lie below the cutoff A* are close to those of H c g. Our result shows that the 
entire operator H restricted to its low-lying energy levels is close to H e g under a slightly stronger 
assumption. 

Theorem A.l Given is a Hamiltonian H such that no eigenvalues of H lie between A_ = 
A* — A/2 and A+ = A* + A/2. Let H = H + V where \\V\\ < A/2. Let there be an effective 
Hamiltonian H eT[ with Spec(ff e fj) C [a,b], a < b. We assume that H e s = n_ff c ffH__. Let D r be a 
disk of radius r in the complex plane centered around z$ = ^±£. Let r be such thatb+e < ZQ + r < A* 

^ 1 A c-^, Hi,; I I'm- rill - CI 



(see Fiaure lXj]) . Let w e g = Kr- Assume that for all z G D r we have — iT e ff|| < e - Then 



uir rr „ ^ 3(||ff cff || + e)||V|| , r(r + z )e 

"<A» — "off S "T 1| rj I, ^7 w N- (A.2) 

A+ - H-Heffll - e {r - w eS )(r - w eS - e) 

Before we prove the theorem, let us make a few comments about how it can be applied. We 
have assumed that H e Q has no support in £ + ; this will be the case in typical applications since 
i7 e ff approximates S_(z) which has support only on It is not hard to modify the theorem if 
H eT[ has (necessarily small) support outside £_ . 

The r.h.s in Eq. (|A.2p contains the energy scale ||-£f c ff|| which is not invariant under shifts by 
aL. In applying the theorem to a Hamiltonian H one can always shift this Hamiltonian H by al, 
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b + e 



A* + A/2 



Fig. A.l. The disk D r in the complex plane, the spectrum of Hes and the other parameters 
in Theorem lA~2l 



without changing its eigenvalues or eigenvectors, such that H e $ has a spectrum centered around 0. 
In that case ||-ff e ff|| = min a ||i? e fr + a A\ = WeS, the effective width. Thus one may replace ||-£^ e ff|| 
by w e g in the application of the Theorem. 

In the construction using mediator qubits, we will choose A_ = and thus A* = A/2. £_ is 
the space in which the mediator qubits are in the state |00 . . . 0) and H e R is of the form H tSLlget (g) 
1 00 . . . 0)(00 . . . 0|. In order for the right-hand-side of Eq. (|A.2|) to be small, we need to take the 
spectral gap A to be sufficiently large (some poly(n)). This will directly bound the first term on 
the right hand side. Now consider the second term and the choice for r. In our applications H e f[ is 
derived from the perturbative expansion of Since H and H on n qubits have norm poly(ra), 

the Hamiltonian H c r (related to the target Hamiltonian) will also have norm poly(n). Hence a, b 
and thus zq are at most poly(n). Note that we need to take zq + r > b + e which implies that 
S_ (z) has to be approximately equal to H e ^ in a range of z which is larger than what is needed 
in Theorem [3 Secondly, it is necessary that the eigenvalues of H are bounded away from A* , the 
difference between the largest eigenvalue below A* and the smallest eigenvalue above A* needs to 
be at least l/poly(n). For our mediator qubit gadgets, one could take r (for example) to scale as 
fal/k £ Qr some constant k > 1 in order for these conditions to be fulfilled. 

Proof, (of Theorem IA.l|) We start from Theorem 3 in Ref. [1] (stated as Theorem lA.il in this 
paper) which shows that under the assumptions in the Theorem one has |Aj(i?<A < ,) — \j(H e s)\ < e 
for each 1 < j < dim(£_). We can draw a contour C in the complex plane, the disk D r in Figure 
lA.lj, that encloses all the eigenvalues of H < \ t and none of the higher eigenvalues of H. The radius 
r needs to be chosen such that b + e < zq + r to include all the eigenvalues of i?<A, • At the same 
time zq + r < A* such that none of the higher eigenvalues of H are included in the contour integral. 
Using Cauchy's contour integral formula we can write 



The remainder of our proof proceeds in two parts. In the first part we show that i?<A» is close 
to n_if < A,n_; this is expressed in Eq. (|A.8|) . In the second part we show that H_ H < \ :t Tl_ is close 
to H e fi, expressed in Eq. (|A.12p . 




(A.3) 
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First part. We have 

||£ <A . - H_ff <v n_ 



= \\u + H <Xt u + + n + H <x n_ + n_if <At n H 
< 2||n + ff <A ,|| + ||^ <A .n + ||, 



(A.4) 



using standard properties of the operator norm ||.||. Let II <Ait be the projector onto the space 
spanned by the eigenvectors of H <Xst with eigenvalues below A*. We can insert tl <Xt before or after 
H <Xt and use that for projectors P±,P2, \\P1P2W = H-P2-P1II) so that 



\\H <Xm - n_# <A „n_|| < 3||# <A J ||n+n <A j| < 3(\\H eS \\ + e ) ||n+n <A j 

In order to bound this, we first derive 

||n+Fn <A j| = ||n + ira+n <A j| > A+||n + n <A j|. 

On the other hand, we have 

||n + #n <A j| < \\u + hu <x j\ + \\v\\ < (||iT eff || + e)||n + n <A j| + 

Putting the last three equations together gives the final bound 

3(\\H e g\\ + e)\\V\\ 



H. 



n_# <Ai n_|| < 



H eS \\ + e) 



Second part. We consider 



2vn 



c 



zh_G(z)n_ dz, 



(A.5) 
(A.6) 
(A.7) 

(A.8) 
(A.9) 



and recall that H-G(z)Tl- = G (z) = (zl- — , Eq. (jA.ip . By showing that this operator 

is close to H-(zI — H c r)~ 1 II- = (zl- — -ff c ff) _1 , we will be able to deduce that H-H <Xt H- is close 
to H ea . 

For all z € D r , — H c g\\ < e by assumption. In order to bound ||(.z/- — S-(z)) -1 — (zl- — 

-f^efr) _1 ||) we will use the following 

\\(A - B)- 1 - A- l \\ = \\{I - A^B)-^- 1 - A- 1 ]] < ((1 - WA-^l \\B\l)- 1 - l) WA' 1 ]], (A.10) 

when \ \B\\ < 1. We choose A = zl- — H e Q and B = — H c r. For z G C (i.e. on the 

contour) 1 1 ^4 1 1 1 < (r — w e fi)~ l and thus ||5|| < T J W < 1. It follows that 

S u V \\(zi--x„(z))- 1 -(zi--H eS )- 1 \\ < e — T . 

zee (r- w eS - e)(r - w eS ) 

Now we will use the following for an operator-valued function F(z) and a contour C with radius 
r around a real- valued zq: 



1 

2ttI 



zF(z)dz 



c 



< r(r + zq) sup ||F(z)||. 

zee 



(A.ll) 



Using this bound and the resolvent for H e ff, we find that 



IL_# <A .n_ - H, 



eff 



< 



J c 
r(r + z )e 



(r - w e g)(r - w eS - e) ' 



(A.12) 
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where we have used that H c q = n_-ff e ffH_. Putting Eqs. (|A.8|) and (|A.12p together gives the 
desired result, Eq. (|A.2|) . □. 

The proof technique used in this theorem can be easily adapted to prove properties of the low- 
lying eigenspace of H ; this is the content of the following Lemma. For the resulting bound of 
Eq. (|A.13|) to be useful one needs (1) to take A large enough compared to ||V|| (which bounds 
the first term in Eq. (|A.13|) ) and (2) the gap A e g- of the effective Hamiltonian H c g, defined as 
A e ff = Ai j6 ff — Ao, e ff) needs to be bounded away from zero. In particular in the Lemma we can take 
r = A ff — 2e and then the second term in Eq. (|A,13j) can be (upper)-bounded by (a ff -2e)(A ff -3e) ' 

If A e ff > po iy( n ) we can thus take a polynomially small e to bound the second term in Eq. (|A.13j) 
by some other inverse polynomial. 

Lemma A.l Given is a Hamiltonian H such that no eigenvalues of H lie between A_ = A* — A/2 
and A+ = A* + A/2. Let the perturbed Hamiltonian H = H + V where V is a small perturbation 
with \\V\\ < A/2. We assume that H e g = n_i? e ffll_ and Spec(H e ^) C [a, b]. Let < e < A and 
assume that for all z G D r , a disk of radius r centered around zq = Ao c fT with e < r < A e g — e, we 
have ||E_(2) — i? e ff|| — e - Let IIo e ff be the projector onto the ground-space of H % g with degeneracy 
d. Let IIi ow be the projector onto the d lowest-lying eigenvectors of H . Then we can bound 

lift n II < 3 W V W i £ ( A Q,eff +r) 

Hlow - IJ-o.cfT S t 7T ; — r H 7 \ — • (A.L6) 

A + - (A ,eff + e) r(r - e) 

Proof. As in the proof of Theorem lA.il we first prove that n_IIi ow n_ is close to ITi ow . Then we 
show that n_IIi ow n_ is close to the projector onto the ground state of H c r, ITo^ff- For both parts we 
will use that due to the assumptions in the Lemma, Theorem [7] implies that for all z = 0, . . . , d— 1, 
|Aj(i?eff) — Aj(-ff)| < e. Let Ao, e fr be the lowest (degenerate) eigenvalue of H c r. We can first bound 

||n low -n_n low n_|| < 3||n+n low ||. (A.i4) 

As before we bound ||n + i2TIi ow || in two different directions: 

( e + A ,eff)||n + ri low || + \\v\\ > ||n+Fn low || > A+||n + n low ||. (A.15) 

These inequalities together with the previous equation give us the first bound 

||n low - n_n low n_|| < - — (A.16) 

In order to prove the other part we draw a circular contour C of radius r centered around 
zo = Ao, e ff with e < r < A e fj — e such that it encloses only the lowest d eigenvalues of H. We will 
choose r such that A e g — r > r or zq + r is closer to Ao, e ff than to Xi >e g ■ We have 

n_n low n_ = -L / n_G(z)n_. (A.17) 

2tti J c 

We use that ||S_(2:) — H c g\ \ < e for z £ C and bound 

sup \\(zl- - ^-(z))- 1 - {zl_ - i^eff)- 1 !! < . 6 . , (A.18) 



zee 

using IK^I- — H c g) < r^ 1 for z G C. It follows that 



r(r — e) ' 



|n_n low n_-no, cfr |l < e(A °' eff + , r) ■ (A.19) 

r[r — e) 



□ . 
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